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We study the Higgs excitation in the Goldstone phase of the relativistic O(N) model in two 
spatial dimensions at zero temperature. The response functions of the order parameter, and its 
magnitude-squared, become universal functions of frequency in the vicinity of the quantum critical 
point described by the Wilson-Fisher fixed point, and we compute them to next-to-leading order in 
l/N. The Higgs particle has an infrared singular decay to gapless Goldstone excitations, and its 
response functions are characterized by a pole in the lower half of the complex frequency plane. The 
pole acquires a non-zero real part only at next-to-leading order in l/N, demonstrating that the Higgs 
excitation has an oscillatory component even in the scaling limit. Both the real and imaginary parts 
of the pole position vanish with the correlation length exponent v upon approaching the critical 
point. We present evidence that the spectral density of the 0(iV)-invariant amplitude-squared of 
the order parameter has a peak at a non-zero frequency in the scaling limit. We connect our results 
to recent experimental studies of the superfluid-insulator quantum phase transition of ultracold 
bosonic atoms in optical lattices. 



I. INTRODUCTION 

Quantum critical points described by relativistic field theories have long been the focus of theoretical study in 
condensed matter physicsP^ A number of experimental realizations of such critical points have also appeared in 
recent years. The critical point of the quantum Ising model in spatial dimension d — 1 has been realized in the 
ferromagnetic insulator^ CoNb206, and also in tilted lattices^ of ultracold bosonsP^ Neutron scattering experiments^ 
on pressure-tuned TICUCI3 have mapped the excitations across the quantum critical point associated with the onset 
of Neel order in a dimerized antiferromagnet in d = 3. The organic insulator (TMTTF^PFg has a quantum phase 
transition^ involving onset of Neel order under pressure in d = 2, but the excitations have not yet been completely 
studied. Ultracold bosonic atoms in optical lattices provide remarkable tunable realizations^^ of the superfluid-Mott 
insulator quantum phase transition in d = 2 and d — 3, whose properties can be studied closely using in situ imaging 
techniques fHHUl 

This paper will study the excitation spectrum in the vicinity of the quantum critical point of the relativistic O(N) 
model in two spatial dimensions. The critical point itself is a conformal field theory in three spacetime dimensions (a 
CFT3), associated with the Wilson-Fisher fixed point Extensive studies of the excitation spectrum have appeared 
in previous workpEHUZI both at T — and T > 0. Our focus here will be on the nature of the amplitude ('Higgs') 
excitation in the Goldstone phase with broken O(N) symmetry at T = 0. It is generally believed that the excitation 



structure of this Goldstone phase is well understood, with N — 1 gapless Goldstone modes which behave like free 
particles at long wavelengths. However, there remain key open questions on the nature of the amplitude mode. In 
mean field theory, there is a sharp, gapped excitation in the longitudinal response, which may be identified with the 
Higgs particle. Beyond the mean field level, this excitation can emit multiple Goldstone excitations even at T = 0, 
and so the stability of the Higgs is not guaranteed. The Higgs has been arguecP^H to be marginally stable in d = 3. 
In d — 2, there are infrared singularities^^ in the spin wave emission process, which indicate that the longitudinal 
response is overdamped, and imply that the ultimate nature of the amplitude mode is a strong-coupling problem. The 
question of whether a finite frequency peak appears at the Higgs energy can depencP^ upon the specific observable 
being used to probe the amplitude mode. 

Here, we will present a 1/N expansion of the spectrum of the Higgs excitation to next-to- leading order, and find 
some qualitatively new features. We show that the Higgs response functions are characterized by a pole in the lower- 
half complex frequency plane, and this pole acquires a non-zero real part at order 1/N. We also find an associated 
structure in the spectral functions on the real frequency axis. Our results can serve as tests of the strong-coupling 
dynamics of the CFT3 in experiments. For instance, recent experiments on ultracold bosons in optical lattices^ have 
studied the low energy response near the 2D Mott-superfluid phase transition, described at integer boson filling by 
an 0(2) critical point with relativistic dynamics in 2+1 dimensions. 

We represent the O(N) order parameter field by the TV-component real field <p a , with a — 1 ... N. Let us orient the 
condensate in the Goldstone phase in the a = N direction, so that (<j) a ) oc S a> jv- As emphasized in a recent paperP 
it is useful to consider two separate observable probes of the Higgs spectrum of fluctuations about this condensate. 

The first observable concerns the longitudinal fluctuation of the order parameter, represented by 4>n, and we denote 
its connected two-point correlator by the susceptibility x, where 

x(p) = J d 3 xe-*p- x [{M*)Mo)) - (M*0XMo)>] (i-i) 

where a; is a spacetime co-ordinate, p is spacetime momentum. We will work largely in imaginary time, and represent 
p = \p\, x = \x\, etc. The real time retarded response function, x( w ) is obtained by the analytic continuation p — > —iu: 
where w is a measurement frequency and the momentum is zero. (Because of the relativistic invariance of the theory 
at T = 0, the susceptibility at non-zero momenta k can be deduced by u> — > \fw 2 — k 2 .) The longitudinal susceptibility 
x( w ) is connected to the neutron scattering cross section for experimental realizations in quantum antiferromagnets. 

The second observable, the scalar susceptibility, is associated with two-point correlations of the O(N) invariant 
amplitude-squared of the order parameter 1 <$? a (in subsequent expression we will omit the summation sign, and 
use an implicit Einstein convention for summation over repeated indices). Specifically, we consider the response 
function 

S(P) = ^J d s xe-^ [<$(*)$(Q)> - {4>l(xWj(0))] (1.2) 

and its analytic continuation S(uj) after p — > —itu. Our main results concern the structure of the function S(uj) in the 
vicinity of the quantum critical point. The recent experiments on ultracold bosons in optical lattices^ measure the 
imaginary part, S"'(w), by determining the energy absorption rate under periodic oscillations in the strength of the 
lattice potential, which induce an oscillation in the hopping matrix element of the bosons. 

The susceptibility x( w ) was computed in Refs. [3HSI In the Goldstone phase, it was found that there was a zero 
frequency peak in the imaginary part x"(w), with x"( w ) ~ 1/w, and no finite u peak which could indicate the energy 
of the Higgs excitation. When considered as an analytic function of a complex frequency Li, the retarded response 
function x( w ) was found to have a pole in the lower-half plane at a purely imaginary frequency, indicating that the 
Higgs mode was purely relaxational. We will compute the leading 1/N corrections to these results in our paper. We 
continue to find that x"(w) ~ 1/u at small u. However, a qualitatively new feature does emerge in the location of 
the pole in the lower-half complex plane: it now appears at a frequency which has a non-zero real part, indicating an 
oscillatory component to the Higgs excitation. 
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In the Goldstone phase, the response function S(uj) also has a pole in the lower-half frequency plane at precisely the 
same location as that for However, now the oscillatory nature of the pole has a more visible signature on the 

real frequency axis. Previous worlP^ found that the imaginary part behaves as S"(w) ~ to 3 at small ui. At larger w, 
no peak-like structure was found in the N — oo computation in the strict scaling limit; this scaling limit is taken by 
including only the single relevant perturbation away from the CFT3. Reference [7] carried out computations beyond 
the scaling limit at N — oo, and found a finite u peak, representing the first signature of a Higgs-like excitation. In 
the present paper, we will compute S(uj) to order \/N. We find evidence for a finite-w peak in S"(u>) even in the 
strict scaling limit. The issue of whether the peak in S"(w) appears in the scaling limit, or not, is important because 
it determines the manner in which the peak frequency, u> p vanishes upon approaching the quantum critical point. A 
peak obtained in the scaling limit will have lo p ~ (g c — g) v , where v the usual correlation length exponent, and g is 
the coupling constant tuning across the quantum critical point at g = g c (here the dynamic critical exponent z = 1). 
Different behavior is expected for a peak which appears only because of corrections to scaling: in this case we can 
estimate from the results in Ref. [7|that w p ~ (g c — g) l '^~' a '', where w is the leading correction to scaling exponent at 
the Wilson-Fisher fixed point P^l Another key difference between the two types of peak is found in the behavior of the 
width of the peak, 7, as one approaches the critical point. For a peak obtained in the scaling limit, the ratio ^/uj p is 
fixed, and hence 7 must necessarily go to zero at the critical point. In contrast, for a peak which appears only due to 
corrections to scaling, 7 will in general not vanish at the critical point. 

While our paper was being completed, we learnt of the numerical study of Pollet and Prokof 'ev.^ They found a 
peak in S"(uj) that appeared to vanish as to p ~ (g c — g) v ; this is further evidence that a peak is present already in 
the scaling limit. 

We will present our scaling limit results for S(ui) (and xC^O) m terms of universal scaling functions. As in the analysis 
of Ref. for the quantum phase transition in the d — 3 coupled-dimer antiferromagnet, it is useful to normalize 
the universal structure in the energy scale by comparing the response functions on systems located symmetrically on 
opposite sides of the critical point. Thus let us take a point in the gapped phase with 0(N) symmetry preserved with 
g = g+ > g c - Then its partner point with g = g- < g c is determined by 

9+ - 9c = g c - 9- (1.3) 

Let A be the exact energy gap to single-particle excitations in the symmetric phase at g — <? + ; consequently there is 
a pole in x( w ) a t the real frequency uj = A for g = g + : 

^ = 1 , 9 = 9 + , (1-4) 

and this defines the value of A. By the definition of the critical exponent u, A ~ (g + — g c ) v upon approaching g c . 
Note that our definition of A is such that the threshold of S"(uj) in the symmetric phase occurs at 2 A; for applications 
to the superfluid-insulator transition, A is the gap in the single boson Green's function in the Mott phase. We will 
use A to set the energy scale of all excitations, including those at g = g—. 

For g = a key signature of the retarded response functions x( w ) an d S(u) is their common lower-half complex 
plane pole. Its location, w po i c , is completely determined by A, and we find 

Wnoic 4 1 /l6(4 + v/21n(3- 2V2)) \ , , , 9k , N 

= + — — i ^ —tL + 2.46531203396 i) +OCL/N 2 ) , g = g- 1.5 

A 7T N \ TT Z I 

The purely imaginary N — 00 result was obtained earlier^. Note the key appearance of a real part ~ 
2.443216943237169/7V in the 1/N correction, which is a measure of the oscillatory component of the Higgs exci- 
tation. 

We have claimed above that x( w ) an d have the same pole structure in the lower-half complex frequency 

plane for g — g_. This is a consequence of the broken O(N) symmetry, which mixes together the fluctuations of the 
longitudinal component and the amplitude-squared of the order parameter. This will be evident from the structure 



3 



of our analysis in Section IV In contrast, in the symmetric gapped phase g = g + , these two observables have very 



different analytic structures. While has a pole on the real frequency axis, as in Eq. (1.4 1, the singularity of S(ui) 
is quite distinct. This response function is in the two-particle sector, and so for g = g + it has a threshold singularity 
at w = 2A, as we noted above. We will show that this threshold singularity is of the form 



In 2 



4A 



2A) 



sgn(w) i 



2A) 



fj = 9+- 



(1.6) 



Here A is a non- universal cut-off dependent prefactor that depends upon the precise definition of <p a , and that can 
be replaced by a constant in the vicinity of the critical point. 

We also compute the universal structure of the amplitude-squared response function, S(u), on the real frequency 
axis in the Goldstone phase with g = g_. Here, there is no threshold, and we have S"(uj) ~ u; 3 as noted above. More 
precisely, we find 



for 



0, g = g-, 



(1.7) 



where A is another non- universal cut-off dependent constant prefactor, and A is defined as always by Eq. (1.4). The 



recent observation of the low frequency spectrum in ultracold atom experiments^ are consistent with the scaling in 



Eq. (1.7). The ratio of the amplitudes of the low frequency singularities at g± in Eqs. (1.6) and (1.7) are universal, 



and we find to order l/N [from Eqs. (3.10) and (4.16 1] that 



0{1/N 2 ). 



For general U), the scaling limit structure of S(co) is (see the end of Sec. [TTJ) 

S(u) ~ constant + A 3 - 2/ "$(w/A), 



(1.8) 



(1.9) 



where $ is a universal function which we will compute here to order l/N; the missing proportionality constant in 
Eq. (1.9) can be fixed by demanding consistency with Eqs. (1.6) and (1.7) for g = g±, and the additive constant is a 



non- universal real number. We will find evidence that in the Goldstone phase (g — g_), the function $ has a peak 
as a function of w, which can be taken as an additional signature of the Higgs mode in two spatial dimensions. For 



small u, the corrections to Eq. (1.7) are in Eq. (4.15) and these are schematically 



S" \uj) ~ io 3 + uj 5 log(w). 



(1.10) 



This non-analytic correction is significant, and a signature of the strong coupling to the gapless modes. Associated 



non-analytic terms are responsible for the real part of the pole of S{u/) in Eq. (1.5), as will become evident from the 
discussion in Section |V] 

We will begin by setting up the general formalism for the l/N expansion for the Wilson- Fisher fixed point in 
Section [nj We will show here that the correlators of 4> 2 a can be efficiently computed by relating them to the correlators 
of an auxiliary scalar field. We will then apply this method to computations in the gapped phase, g = g + , in 



Section III Our main results for the amplitude excitations of the Goldstone phase, g 
present a unified review of our results in Section W\ 



<7_ , appear in Section IV We 
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II. GENERAL FORMALISM 



In this section, we set up some general formalism for the O(N) model. We will be interested in the field theory 



Z U [J] = J Vcj) a exp (-J 



^f + ^^-N/gf + Jcg 



(2.1) 



where J = J d 3 x, u is a quartic non-linearity, and g is the tuning parameter across the transition. We have included 
a spacetime dependent source for the amplitude-squared mode, J(x), and will eventually set J(x) = after taking 
functional derivatives of Z U (J). Note that the partition function in Eq. (2.1) is just the standard </> 4 field theory of 



the Wilson-Fisher fixed point, in an unconventional notation; the notation is designed to simplify the large N limit, 
and to make a natural connection to the fixed length non-linear sigma model in the limit u — > oo. 
We decouple the quartic term, and write the partition function as 



Z U [J] = / Vcj) a VXexp - 



(5</> Q ) 2 + JLA (4>l ~ N/g) + ^ + 2 .m 



(2.2) 



In the large N limit, we perform the Gaussian functional integral over <f> a (at J = 0), and expand the fluctuations in 
A about the saddle point of the effective action. In the symmetric phase, the saddle point is at iX/y/N = r, where 



1 

- + 



2u 



1 



(2.3) 



where / = f d 3 p/(8ir 3 ). The nature of the saddle point in the Goldstone phase will be discussed later. So we write 
A = —iyNr + i2vN J + A, and completely integrate out the <fr a field. Then, the effective action for A is 



Z U [J] 
S 



VXexp (~Sj — Sq — Si) 



S.j = 



i_ r /%) 

Nr Nr 2 



l 

4u 



8 a 



J 



Nr + 2Nu/g + i\/NX\ NJ 



2u 



2u 



(2.4) 



where Si involves cubic and higher order terms in A alone, and is specified in Appendix [Sj we note that Si is 
independent of u. We show in Appendix [B] that r plays the role of the tuning parameter away from the quantum 
critical point, with r ~ (g — g c ) 2 for g > g c . 

Now we can take functional derivatives with respect to J, and obtain the two-point correlator of § 2 a at J = as 



£(*)$(0)) - (cbl) 2 = *5 d (x) - [(A(x)A(O)) - (A) 5 



(2.5) 



So we see that the amplitude-squared correlator, S(p), is simply proportional to the two-point correlator of the A 
field. Specifically, we will identify 



S(w) = constant — G\\(ui), 



(2.6) 



where the right hand side contains the connected A Green's function, and we are not interested in the additive non- 
universal constant. The function G\\(uj) will be computed in this paper in the 1/N expansion; this expansion can be 
carried out at any value of u, including the theory in which we impose the fixed length limit by taking the u — > oo 



limit. We also note that the identity (2.5) also holds in the Goldstone phase. Having proved Eq. (2.5), we will set 
J = is all subsequent computations. 
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In the symmetric phase, g > g c , we can write the A Green's function as 



[GAACrir^IGM^-ExxCp). 



(2.7) 



The bare propagator is 



where 



IL(p,r) + l/(2u)' 



(2. 



n(p,r) 



tan 



2^ 



(<? 2 +r)((p + q) 2 + r) 4irp 

= -^lnf 2 f + tP 
8717W \2vr — ip 



(2.9) 



Note that r) diverges as p,r -> 0, and so the l/(2u) term in Eq. (2.8) can be neglected in this limit. Because 
u does not appear in the non- linear terms in Si, we need to send u — > oo to obtain the scaling limit, and we will do 
so in most of our expressions; thus the scaling limit response function in the vicinity of the quantum critical point 
is the same for the "soft spin" (finite u) and fixed length (u — oo) theories. In terms of the expressions obtained in 
Appendix |A| the expression for the self-energy to order 1/N is 



Eaa(p) = J [K 3 (p,k,\p + k\)} 2 G xx (\k + P \)G xx (k) 



1 

2iV 

1 
6N 



K 3 (p,- Pl 0)G° xx (Q) / K 3 (k,-k,0)G° xx (k) 



^ f Ki {p,k-p,-k)G° xx {k) + ^ f K A {p-p,k,-k)Gl x 



(k). 



(2.10) 



We can also obtain the expression for the 4> a Green's function, G(p), by similar methods. To order 1/N we obtain 

1 



G-\p) = P 2 + r+^J G° xx (q) 



1 

2N 



G° xx (0) 



((p + q) 2 + r) 
1 



k, q (k 2 +ry((k + qy + r) 



G° xx (q)- 



(2.11) 



We describe the structure of the critical singularities, as g approaches g c , in the expressions in Eqs. (2.101 and (2.11 ) 
in Appendix |B| 

Before proceeding further, we give a brief derivation of Eq. ( 1.9 ), which plays a central role in our analysis. Consider 
the D — 3 dimensional statistical mechanics model defined by the partition function in Eq. (2.1 ), and define the free- 
energy density / = — y In Z U [J = 0]. Then, taking two derivatives of / with respect to g _1 we obtain, 

Q2 f 



= -uN - 



u 

V 



(2.12) 



which is the p = Matsubara frequency correlation function S(p = 0), up to additive and multiplicative constants. 
On the other hand, near the critical point g Cl / can be split into regular and singular parts, / = f reg + fsing where, 
according to the hyperscaling hypothesis, f S i ng ~ £,~ D - Here, the correlation length £ ~ \g — g c \~ u or, equivalently, 
i ~ I <?cT 1 — 9 l \ " • -By taking two derivatives of f S in g with respect to g~ x , we conclude 



S(p = 0) - £- D + 2 /» + (regular constant). 



(2.13) 



G 



Finally, since £ oc 1/A for a relativistic theory, and generalizing this result to arbitrary real frequencies, we obtain 



S(w) ~ A 3 ~ 2/l/ $(w/A) + constant, 



(2.14) 



which is Eq. ( 1.9 ) 



III. SPECTRAL FUNCTIONS IN THE SYMMETRIC PHASE 



This section limits attention to the gapped phase where the O(N) symmetry is preserved, present for g > g c . The 
dominant singularity in the </> Q propagator is the pole corresponding to the TV-fold degenerate particle excitation. Let 
us determine the position and residue of this pole to order 1 /N. 



The quasiparticle pole is determined by the solution of G (p) = where from ( B9 ) 



G-\p) = p 2 



dq 



NttpJq ^tan- 1 ( (? /(2^)) 



In 



(p + q) 2 +r\ 4pq 
(p — q) 2 + r J q 2 + 4r 



(3.1) 



So the pole is at p 2 — — A 2 where 



a 2 = r H q2 

l+ N7rJ aq t^\q/{2^r)) 



tan 



The integral yields 



A = r 



32 / f A 2 
m 



3vr 2 Ar 



0.50260886161586 



(3.2) 



(3.3) 



This is the needed connection between the renormalized energy gap, A, and the bare tuning parameter r. Recall that 
we showed in Appendix LB] that r ~ (g — g c ) 2 , and so A ~ (g — g c Y , where 



v = 1 



32 
3tt 2 N' 



Also, we have^ for the wave-function renormalization Z: 



1 

Z 



dG-Hp) 



dp 2 



1 



A 



dq- 



= 1 + r) In 



Nnr J y tan- 1 (q/(2Vf)) 
A 



' tan- f 2 ^ 
s/r \ q 



Vr 



1.78297175565687) 



q 2 + 4r q 



(3.4) 



(3.5) 



where the critical exponent r\ is given by Eq. ( B6 1 



A. Threshold singularity 



As we discussed in Section [TJ the low frequency singularity in G\\ is a threshold at a frequency 2 A. So we consider 
the non-analytic terms in G\\(p) at p — — i(2A + lo) as uj — > 0. 
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At N — oo we have 



Tl(-2iy/r- 



In 



4v^ 



iw, r) 



(3.6) 



So near threshold we have 



ImG° AA (p) = 



32w 2 ^r 



ln 2 (4 v ^/|u;| 



(3.7) 



The 1/ln 2 prefactor is easy to understand. We are considering the creation of two slowly moving particles in two 
dimensions with a k 2 dispersion at low momenta; such particles have a well-known logarithmic singularity in their 
T-matrix for short-range interactions.^ This argument also implies that the 1/ln 2 prefactor will be robust at higher 
orders at 1/iV, as we indeed find below. 

At order 1/N, we compute the threshold singularity in G AA from Eq. (2.10|. Notice that the second and fourth 
terms in Eq. (2.10) can be absorbed into a self-energy renormalization of the <p a propagator, which only contributes 
via the quasiparticle wavefunction renormalization to the threshold singularity. So for the purposes of the threshold 
singularity we can write 



G\x(p) 



[K 3 {p,k,\p + k\ 



Z 2 n(p,A 2 ) NTP{p,r) J k IL(\k+p\,r)Il(k,r) 3NIl 2 (p,r) J k 



K 4 (p,k,-p,-k) 
II(fc,r) 



(3.8) 



It remains to evaluate the k integrals in the limit lu — > 0. We present details of this evaluation in Appendix [C] The 
final result can be written as 



AA3-2/u 

ImG AA (p) = 6(u) 
In (4A/|w|) 



(3.9) 



with the amplitude A is that appearing in Eq. (1.6), and is given by 



A = 32tt 2 A 2 ^- 2 1 



1.481 

N 



(3.10) 



IV. SPECTRAL FUNCTIONS IN THE GOLDSTONE PHASE 



We write <p a — (vNao + a, 7i"i) with i = 1 . . . N — 1, where 

o i r i 



~ ~ / 2 • 

9 JpP 2 



(4.1) 



is the condensate at N = oo. From Eqs. (1.3) and (B7), we can relate this to the value of r at the partner point in 
the symmetric phase 



4-7T 



(4.2) 
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where r is related to the single particle gap in the paramagnet via Eq. (3.3 1. This yields a relation between a 2 and 
the gap at the partner point in the symmetric phase: 



324 ln 



u 4tt 3tt 2 N 
where C = 0.5326726500732776. This can also be written as 

C 



A 
16^ 



Col 
N 



A 

4tt 



N 3n 2 N n ^4A J 



within the same order in 1/N. 

We can set A = A in Eq. (2.2), and then, integrating over the 7Tj we obtain the partition function 



Z = J VaVXexp (S -S x - S 2 ) , 



So = 



S 2 = — 



p 2 a 2 + 2icr a\ + -U(p, 0)A 2 



1 



2VN \9 



A- 



4N 



n( P ,o)A 2 



2VN 



Act 2 , 



(4.3) 



(4.4) 



(4.5) 



where Si is just as in Appendix [X] but with K 3 4 evaluated at r = 0. Now the connected a-a Green's function is 
equivalent to the longitudinal susceptibility x( w ) introduced in Sec.|l] 



(4.6) 



The bare connected Green's functions of Eq. (4.5) are 



gUp) 
gIAp) 



n(p,o) 



p 2 II(p,0) + 2a 2: 
2p 2 

p 2 U(p,Q) + 2a 2: 

—2iuQ 
p 2 U(p,0)+2a 2 - 



(4.7) 



The corresponding renormalized connected Green's functions require computations of the self-energies at order 1/N, 
and these are presented in Appendix |Dj 



Using the value of H(p,0) in Eq. (2.9), we see that all the Green's functions have poles at p = — 16<7q, which 



corresponds to the lower-half plane value u = —HQa 2 ,. After including the self-energies, this pole corresponds to the 
point where the effective quadratic action for A and a has a zero eigenvalue. Equivalently, this is the condition that 
the determinant of the quadratic form vanish, i.e., 



det 



P 



^aa(p) 



iljQ - Y.<j\{p) 



ia - ^a\(p) l/(16p) - Saa(p) 



(4.8) 



where S are self-energies computed in Appendix D Note that upon ignoring the self-energies, the solution of Eq. (4.8 ) 



is indeed p = — 16a 2 . Rather than solving this equation with self-energies to next order in 1/N, it is more convenient 



to locate the pole by examining the zero of 1/G aa , which we will do in Sec. IV B 



A. Computation of $ 



We can now use the results for the self energies computed in Appendix [D] to obtain expressions for the connected 
A Green's functions to order 1/N: 

Gxx( P ) = G° xx (p) + Tj\\(p) [G° AA (p)] 2 + E„(p) [GUp)] 2 + 2^(p)G xx (j>)G° X(r (j } ) 



I6p 2 



I6p 3 



1024 p 2 (p + 8ct 2 ) 



In 



A 



p+lQal A(p+16cr 2 ) 2 3tt 2 N( P + 16a 2 ) 2 \16a 2 
256cr 2 



(p + 16a 2 ) 2 N 
256A 8192cr 2 / A 



7T 2 A 



TT 2 N 



[p 2 Fxx(p) - 4F°M + 2pa 2 F aX ( P )] 
In 



V 16a 2 ' ' 



(4.9) 



where the F functions are defined in Eqs. (D2 |-|D4|). The terms in the last line are p-independent real constants, and 
will be dropped since they do not contribute to S"(oj) upon using Eq. (2.6 1. After this subtraction, it is notable that 
all the terms linear in A have canceled out of this expression, indicating that the S(uj) is universal near the quantum 
critical point. 

In order to extract the universal scaling function $ from G xx , as introduced in Eq. (1.9), we first collect the p- 
independent terms appearing in the last line of Eq. (4.9 1. These terms are explicitly real and can be absorbed into 

1024p 2 



the constant in Eq. (1.9). We then use Eq. (4.4) to substitute ao by A, and find: 

16p 2 (p + 4CA/7r) 



G\\(p) = (const) + 



16p 2 



In 



7rA 



P + AA/tt N( P + 4A/tt) 2 3tt 2 N {p + 4A/tt) \AA 
4A r 16tt 2 p 2 F xx (p) - A 2 F aa (p) + 8npAF aX (p)] 



Nir 3 (p + 4A/ir) 2 



Using 



this can be written as 



where 



(A/A) 



2-2/v _ 



= 1 



64 
3tt 2 A 



MA/A), 



G XX ( P ) = (const) + A (A/A) 



2-2/v 



Hx( P /A) 



I6p 2 



l6p 2 (p + 4CA/ir) 



Hx( P /A), 



1024p 2 



(4.10) 



(4.11) 



(4.12) 



/tt) V4/ 



A(p + 4A/tt) AA(p + 4A/tt) 2 3tt 2 N (p + 4A/tt) 

4 [Wn 2 p 2 F xx (p) - A 2 F aa (p) + 8wpAF aX (p)] 



(4.13) 



A7r 3 (p + 4A/7r) 2 

Importantly, notice that all the logarithmic dependence upon the cutoff A has dropped out, yielding a universal 
expression for H x ; this is a significant check on our computation. By comparing Eq. (4.12) to Eqs. (1.9) and (2.6), 
we see that the scaling function $ is 



= -H x 



— lUJ 



(4.14) 



Hence, $ is obtained from the analytic continuation of H x to p 
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FIG. 1: Universal component of the scalar spectral function, 4>"(cj/A), for different values of N. The change of sign in the 
N = 2 curve at large uj/A is an artifact of the approximation being used. 



At small values of z = uj/A, the expansion (D5| yields 
$(z) = 



25.898 \ 2 

47T ] Z + I 

N J 



ITT 



50.584i 
N 



7T 
2lTT 2 

~3N 



25.842 

N 



7T 

T 



40.690 \ 4 



N 



]n(-iz)j z 5 + (z 6 ) . 



(4.15) 



Note that the imaginary part of <!> has a leading uj 3 behavior. This relies on a precise cancellation of a term ~ uj, and 
is consistent with previous results^!. The coefficient of uj 3 is a measure of the decay rate of the Higgs excitation into 



two gapless Goldstone modes. Comparing the definition in Eq. (1.7) with Eq. (4.15) we have 



A = K 2/v - 2 



7T 2 - 



25.842 
N 



(4.16) 



Note that there is a branch cut in Eq. (4.15) at u> = Q, as indicated to lowest order in uj by the term uj logo;. This is 



a non-analytic consequence of the presence of gapless excitations. This same non-analytic term is responsible for the 
real part of the pole in Eq. (1.5), as is shown in Appendix [P| and discussed in Section [Vj 

We next evaluate $ numerically at general values of uj. For this, we must perform an analytic continuation of the 
integrals F in Eqs. ( D2|D4 ) to real frequency p = —iui. A Wick rotation of these integrals can be performed by taking 
p = e~ la ui and k — e~ ia v, and varying a continuously from a — (the Euclidean time result) to a = tt/2 (real time). 
On general grounds, we expect not to cross any singularities upon this Wick rotation — indeed, this can be verified 



explicitly by inspection of the integrands. Thus we interpret F(p = —iui) in Eqs. (D2)-(D4) as integrals along straight 
line segments on the complex plane, running from k — to k = —iui, and from k = —iui to k = —ioo. 

The universal scaling contribution to the scalar spectral function, $"(a;/A), is shown in Fig. [l]for N = 2, N = 3, 
and N — 4. Note that <J>" has a peak in all cases. This is one of the main results of our analysis. It supports 
the appearance of the Higgs mode as a finite width peak in scalar measurements in the quantum critical region. 
The presence of this peak is substantiated by recent experiments^! and quantum Monte-Carlo simulations^ of the 
Mott-superfluid transition at integer filling of bosons, corresponding to the case N = 2. 

The N — 2 curve in Fig. [l] is negative at large uj, and also at very low values of uj, as can be inferred from the 
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FIG. 2: 



Universal component of the scalar spectral function $(cj/A) for N 

For a 



Eq. (4.181 for different values of a 
U /A. This unphysical result is absent at a = 0.4 and a = 0.5 



2, using for N = 2, using the Pade approximant 
0, corresponding to Eq. (4.171, $" becomes negative at very small and very large 



sign of A in Eq. (4.16) for N — 2 (resulting in a very shallow negative region that cannot be discerned in the scale of 



Fig. [I]). These regions of negative <&" are unphysical, since the spectral function must always be positive. Note that 
the scaling function can be written as 



(4.17) 



The change of sign in $" occurs because the 1/N correction $i can become dominant over the N — oo result, $0j 
both at very low and very large z. This is an indication of the breakdown in the 1/N expansion at the current order 
in the case N = 2. One possible method to avoid the unphysical regions is to introduce a Pade approximant, 



Pade 



» 



$ (Z) + ^$!(Z) 

$i(z)/$o(*) 



)(l/iV 2 ) 



(4.18) 



N 



which is equivalent to Eq. (4.17) to our working order in 1/N. Then, by choosing an appropriate value of a 0, 
one obtains a positive spectral function, as shown in Fig. [2} For instance, in the small z region, Eq. (4.16) would be 
replaced by, 



A 



Pade 



k 2/u- 



r 2 - 25.842(1 -a)/N 
1 + 25.842a/(7r 2 7V) ' 



(4.19) 



which matches Eq. (4.16) to our working order in 1/N, yet is positive for all values N > 2 provided a > 0.2366. 



The breakdown in the 1/N result in fact occurs for all iV: at large enough z, <I>" becomes negative for any N, 
although the value of z where this occurs increases with N. This means that there is significant ambiguity in the 
large-z dependence of as seen in Fig. [2| However, many features of $" are robust. For example, $" displays a 
peak provided that the value of a is moderate. For a = 0, the peaks in Fig. [l] are located at w poa k f=s 2.5A for N = 2, 
Wpeak ~ 3.5A for N = 3, and w p0 ak ~ 4.5A for N = 4. There is uncertainty in the locations of the peaks, since these 
do shift when a is varied. However, the fact that the location of the peak in $" is monotonically increasing with N 
is robust for any value of a. 
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B. Computation of longitudinal susceptibility 



The longitudinal susceptibility is related to G aa via Eq. (4.6). From scaling arguments we expect this response 
function to satisfy 



G-» = A 2 -"$ CT ( W /A). 

We will next compute the universal function $ CT to order 1/N. 
To order 1/N, 



G aa {p) GUP) 



P (p+l6a 2 Q ) 



16pa 2 t 8p(p + 80cr 2 ) f A 



AT 



3tt 2 AT 



V 16a n 2 



^ - 256F AA (p) + 32F CTA (p)] , 

where the functions F are defined in Eqs. ( [D2p4| . Using Eq. Q, we find: 

1 . „ ... 4(1 - C)pA 8 P (p + 4A/tt) , 

p (p + 4A/vr) + _J F + FKF _ 7A ' log 



ttN 



3tt 2 N 



7rA 

4A 



A 2 



16n 2 N 



[F a<T (p) - 256F AA (p) + 32F ffA (p)] 



A 2 (A/A)-"tf ff ( P /A), 



where we have used Eq. ( B6 ) to write 
and where we have introduced 



(A/A)-" = 1 - 



8 



3tt 2 AT 



ln(A/A) 



p(p + 4A/ir) 4(1 -C)p 8p(p + 4A/7r) / tt 

A~2 ' _ AT A ' o_9 A 9 AT 



1 

16tt 2 AT 



ttNA 3tt 2 A 2 A' 
[^(p) - 256F AA (p) + 32F ffA (p)] 



which is seen to be a universal function of p/A. Comparison with Eq. (4.20) gives 



At small values of z = w/A, we can use Eq. (D5) to obtain 

_ , , /4 0.5119N . 

0.2793i 2Ai 



1 



2.2811 \ 



A r 



97rAr 



N J 
ln(-iz) ] z 3 + C(z 4 ) 



(4.20) 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



Note the linear dependence at low z, corresponding to the infrared divergence in the longitudinal susceptibility 
X"(w) ~ 1/wP™ 

We compute $ ff numerically performing an analytic continuation as before. Figure [3] shows the imaginary part 
of l/$cr for the cases N = 2, A" = 3, and AT = 4. This corresponds to the universal scaling form of x"(w). In all 
cases, the curves are dominated by the infrared divergence ~ 1/u). In order to suppress this infrared divergence and 
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FIG. 3: The imaginary part of l/<E> CT (o;/A), obtained from the analytic continuation of Eq. ( |4.24| | for N = 2, N = 3, and N = 4. 



This corresponds to the universal scaling form of the longitudinal spectral function x"(^)- Inset: uj/A ■ Im(l/$ a 
TV = 3, and N = 4. 



for N = 2, 



highlight the features appearing at intermediate values of w/A, we plot uj/ A ■ Im(l/<fr a ) in the inset of Fig. [3] . 

Equations (4.22| and ( 4.24[ ) display a key feature of the longitudinal Green's function, which first appears at our 
current working order, 0(1/N). For N — oo, the longitudinal Green's function has poles at p = and p = —AA/tt. 
The latter pole, which lies on the negative p axis, corresponds to an overdamped response with purely imaginary 
frequency uj = —MAjix. Remarkably, at 0(1/N), this pole acquires a real component. The appearance of a real part 
in uj leads to a real-time longitudinal response which has both decaying and oscillatory components. 



In order to find the shift in the location of this pole, we need to solve Eq. (4.8) to order 1/N. Equivalently, we can 
look for solutions to H a {p — — 4A/-7T + Sp) = 0, where Sp = 0(1/N). This gives 



Sp = 



4A 



1-C 



1 

256 



(F aa + 32F aX - 256F AA ) \ + 0(1/N 2 ) 



(4.27) 



where the functions F are to be evaluated at p = —AA/ix. This is carried out in Eq. (D7|, where it is found 
that F\\{— A/w) has an imaginary component - this is the source of the real frequency shift in the pole, u) po \ c = 
i(-4A/TT + Sp). We find 




2.46531203396 i 



0{1/N 2 ) 



(4.28) 



The real part is 16 (4 + \/21og (3 - 2\/2)) /(/Vtt 2 ) - 2.443216943237169//V. This pole, which is most easily extracted 
by studying the longitudinal susceptibility, will also appear in the scalar response function $>(z), as discussed below. 



CONCLUSIONS 



Our paper has studied the nature of Higgs excitations of the relativistic O(iV) model in 2+1 dimension in the 
Goldstone phase in the vicinity of the Wilson-Fisher fixed point P3 It is remarkable that after nearly 40 years of study, 
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zero temperature characteristics of this venerable fixed point have still remained uncovered. 

Our main results were for the spectral function of the amplitude-squared of the order parameter that is measured 
in recent experiments on the supcrfluid-insulator quantum phase transition in Ref. I19[ and in numerical simulations 
in Ref. 25 We expressed our results in terms of the universal complex scaling function $ in Eq. (1.9), which measures 



frequency in terms of a single energy scale A which vanishes at the critical point with exponent v. We obtained 



information on distinct features of this function: it was expanded in powers of 1/N in Section IV A expanded at low 



frequencies in Eq. (4.15), and we determined the location of the pole in the lower-half complex frequency plane in 



Section IV B The pole is specified in Eq. 1.5 note that both its real and imaginary parts are proportional to A, and 
so the width of the Higgs particle is of the same order as its energy. 

Here we combine these results to obtain useful interpolating forms for comparing with numerical or experimental 
data. We consider $ as function of the complex variable z — ui/A — ip/A. An explicit form with all the necessary 
analytic features is 



1 + ic\Z + C2Z 2 + ic^z 3 + ic^z 3 ln(— iz) 



Co 



(5.1) 



where cq. 1,2,3,4 are real constants. The constraints on these constants are: 

• $(z) is real on the imaginary axis in the upper-half plane, z — ip/A with p > 0. 

• §{z) is analytic in the upper- half plane. 

• There is a branch cut in $(z) starting from z — 0, we take this cut to lie in the lower-half plane. Upon analytically 
continuing <&(z) clockwise from the imaginary axis in the upper-half plane we find a single pole in the lower-half 



plane. The position of this pole should coincide with Eq. (1.5). On the other hand, if we analytically continue 



$(2;) anticlockwise from the upper-half plane we again find only a single pole in the lower- half plane, but its 
real part has changed sign. One convenient choice is to place the branch cut along the lower imaginary axis, 
leading to the analytic structure shown in Fig. [4] This analytic structure is remarkably similar to that found in 
recent studies using the AdS/CFT correspondence!^!!!! These studies describe a particle interacting strongly 
with a gapless continuum of excitations, and so bears a similarity to our study of the Higgs particle interacting 
strongly with the Goldstonc excitations. 



• The small z expansion of $(z) should agree with Eq. (4.15). 

v 3-2/v j n ^jj_ e ii m it z» 1 5 as required in order to ensure that S(uj) is independent 



z 



In addition, we expect $(2) 

of A when uj 3> A. Equation (5.1) does not satisfy this requirement. However, a comparison with experiments or 



with numerical simulations is bound to be complicated at high frequencies by lattice-scale contributions. Thus, we 



choose to focus on the regime in which the frequency is not much greater than A. Then Eq. (5.1 ) can serve as a useful 
prototypical form to fit the data. 



Finally, we note that the scaling function ^ a (z), which specifies the response function x(w) via Eq. (4.20), has a 
very similar analytic structure. The only change from $(z) is that there is an additional pole at z = 0. However, the 
locations of the poles in the lower-half plane are the same, and there is also a branch cut starting at z = 0. 
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FIG. 4: Analytic structure of the retarded Green's function $>(z) in Eq. (5.11. This is analytic and unique in the upper-half 
plane, but its structure in the lower-half plane depends upon the choice of a branch cut originating at z = 0. If we analytically 
continue from the upper-half plane, we find a pole in the lower-half plane, but its location depends upon the manner in which 
we encircle z — 0. The figure shows the choice of the branch cut along the negative imaginary axis. Then the pole is at 
Eq. ( |4.28[ ) or at the partner point with the opposite sign for the real part. Note the similarity to the analytic structure of the 



strongly coupled gapless systems studied using the AdS/CFT correspondence in Refs. [27H29 



Appendix A: Nonlinear terms in action 



This appendix will specify the non-linear terms, <Si, in the functional integral over A in Eq. (2.4). They are given 

by 

Si = / K 3 (pi,p 2 ,p 3 )\(pl)\(P2)KP3) S 

"V-1V Jpi,P2,P3 \i=l ) 

l^f #4(pi,P2,P3,P4)A(pi)A(p 2 )A(p 3 )A(p4)<5 ( ^2 Pi J • (Al) 



24N 

The three-point vertex is given bjESl 

#3(Pl,P2,P3) 



1 



(q 2 + r)((q + piY + r)((q - p 2 ) 2 + r) 



1 -1 p 

tan 



AttP ^/r{jp\ + pi + pi + 8r) 

= K 3 { Pl ,p2,p 3 ), (A2) 

where 

P 2 = vlvlvl + r(p 2 i +p 2 2 + pj) 2 - 2r(pi +p\+ p\). (A3) 
At the critical point, with r — > 0, we have 

K 3 { Pl ,p2,p 3 ) = ^— ■ (A4) 
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At the critical point, we will also need 

11(0,0) f 1 ( 1 1 



<z 4 V(g + p) 2 p 2 



n(o,o) 

n2 



where we have used the curious fact that 

1/1 1 



IpP 4 \(p + q) 2 q 2 

only in 3 spacetime dimensions. Similarly, the 4-point vertex is 



Ki(p,k, -p, -k) = 3 / 

Jq 



(q 2 + r)({q + p) 2 + r)((q+p + k) 2 + r)((q + k) 2 + r) 

q 2 +r + (q+p + k) 2 +r - (q + p) 2 -r - (q + k) 2 - r 



and from Ref. |3"T1 at the critical point, 

1 



K 4 (p,-p,k,-k) = 3 / 

J Q 



q 4 (q + p) 2 (q + k] 



(A5) 



= (A6) 



K±{pi,Pl,P3,Pi) = 3 / ■ yz- — -7- — — — ■ r^— — r^— — -. (A7) 

Jq (l +r)((q + Pi) 2 +r)((q + p 1 + p 2 ) 2 + r)((q - p 4 ) 2 + r) 

We will not need this expression in general, but only 2 special cases; we have 

1 



2(p ' k) J q (q 2 + r)((q + p) 2 + r)({q + p + k) 2 + r)((q + k) 2 + r) 

T \- [K 3 (k lP , \k-p\) - K 3 (k,p, \k + p\)} , (A8) 
(p-k) 



p 2 k 2 J q q 4 \(q + p) 2 {q + k) 2 p 2 k 2 

3 M + 3(p • fc ) (A9) 

p 2 k 2 Sp 3 k 3 \p-k\' { ' 



Appendix B: Critical singularities 



This appendix will review the derivation of the quantum critical singularities at g = g c using the basic expressions 
in Section |IIj We will set u = oo in all expressions here. 

Let us first determine the value of the critical coupling g c to order X/N. The critical point is determined by the 



condition G 1 (0) = 0. So at the critical point r = r c where, from Eq. (2.11), 

2 r l l 2 i r l 

Tc = ~Nj q U(q, 0) q* + N n(0, 0) 7 piq P 4 (p + q) 2 U(q, 0) 

2i r i r i / i i 



i 



iVn(0,0) J q Tl(q,0)J pP 4 \(p + q 



(Bl) 



Note that 11(0, 0) is infrared divergent, but we assume this is suitably regulated by putting the theory in finite box; 

. Using Eq. ( A6 ) , we conclude that 

r c = 0(1/N 2 ). (B2) 



we will find that 11(0,0) will eventually cancel out. Using Eq. (A6), we conclude that 



17 



So from Eq. (2.3|, we obtain 



1 



1 



0(1/N 2 ). 



9c J p T 

Next, we determine the exponent r\. For this we need G~ 1 (p) at the critical point r = r c . This is given by 

2 f 1 / 1 1 



(B3) 



n(g,0) \(p + q) 2 q 2 



We have 



where 



G- L (p))=p 2 1 + 7? In 



'7 



A 



3tt 2 A/" 

For the exponent, f, assume the coupling is g, and define r s by 



1 1 



1 1 



9c g J p \P 2 P 2 + r g 



Then from Eq. (2.3) and (B3) we have 



4tt 



r = r g + 0(l/N 2 ). 



Now combining Eq. (2.11| and (B8), we have the result for the inverse susceptibility 

G~ l (p) =p 2 +r + F(p,r) 

where 



F(p,r) 



1 



NL n(«,r) 



1 



%) 

((p + <z) 2 +r) 2n(0,r) 



with 



!T(p t r) S aiI ^ r ^ -2 



Or J q (q 2 + r) 2 ((p + q) 2 + r) 47ryV(p 2 + 4r) ' 



A useful check is to note that G 1 (0) = as r — ► 0. From (BIO) we have 

rA „3 



F(0,r) 



12r 



Mr Jo (g 2 + r)(g 2 +4r)tan- 1 (g/(2Vr 7 )) 



From this we obtain the exponent 



12 /A 2 



7 = 2 



24 
Ntt 2 ' 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 
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Let us now obtain the connected Green's function, G\\(p) to order 1/N at the critical point r — r c . We insert 
( |A4|A5|A8|A9| into ( |2.8|2.10[ ). All occurrences of 11(0,0) cancel out, and 



G xx (p) = I6p 



512 



N J k k\k-p\ 



256p 
N 



(p-k) 



\k-p\ k+p\ 



512 r (p-k) 

~W j k pk 2 \ P -k\ 



128 . 256n /A 



256p /, /A 
In 



256A 

' Nn 2 



16p 1 



P 

64 /AN 22 
3tt 2 N n \p) + Y 



3Nir 2 



The scaling dimension of A is the same as that of </> 2 , which is 3 — l/v, and so we expect 

G\\(p) ~ constant + p 3 ~ 2/u 



and so (B14) is consistent with (3.4) 



(B14) 



(B15) 



Appendix C: Computation of threshold singularity 



We describe the evaluation of the k integral in Eq. (3.8 1 at p = —i(2A +a>) as to — > 



The singularities arising from the integral of fc in Eq. (3.8) are all at energy 4A, and so we only need the explicit 



singularities of all the terms as p — > —i(2y/r + uj). We list these: 



II(p,r) 
n(|p + fe|,r) 
K 3 (p,k,\p + k\) 
Ki(p,k,-p,-k) 



1 



In 



16ir^/r \ —lu 

l + (l + ikx/^F-k 2 /(4r)) 1 / 2 
16^^(1 +ikx/y/r- k 2 /(4r)y/ 2 V 1 - (1 + ikx/y/f - fc 2 /(4r)) 1 / 2 



In 



1 



167r^/rfc(fc — li^frx) 
3 



In 



In 



167rrfc 2 : 



Tm 



k^Jr(k — 2iy/rx) 2 
(fc 2 + 4r)(fc - Aiy/rx) 
1 



In 



(fc - 2i«/rx + (k 2 - Qik^/rx - 4r(l + x 2 ))w/(2Vr(fc - 2iVrx))) 
4 \/''\ /" k^/r(k - 2iy/rx) 2 



In 



(fc 2 + 4r)(fc - 4iy/rx) 



(CI) 



where x = (p-k)/ (pk) 



For the second term in Eq. (3.8), we can now pick out the co-efficient of the threshold singularity to be 

k^/r(k — 2iy/rx) 2 



n 2 N 



In 



- 1 ! " x 



1 



32lTy/f 

N 



hi 



4^ 



doc I dk In 

2 7_i Jo V(^ 2 + 4r )( fc - 4*^)7 (fc- 2iVr ; x) 2 II(|fc + p|,r)n(£;,r) 

l 

x 4.507333003 



(C2) 



For the last term in Eq. (3.8), the term associated with the explicit hx(A^/r/(— u>)) term in K4 is 

128fc 



1 

N 



In 



/4y^ 



V — w 

x Im 



, dx I dk— , 
2./_i Jo xtan- 1 (fc/(2V^)) 

1 



(fc - 2iVfx + (fc 2 - Qiky/rx - 4r(l + x 2 ))uj / (2y/r(k - 2iy/rx))) 



32ny^ 
N 



In 



4^ 



2 In 



4^ 



-A-r 



3.44906254 



(C3) 
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The ln(A) term above has the same co-efficient as that expected from the corresponding term in Eq. (Bf4l 



The remaining contribution of the last term in Eq. (3.8) is 



f 

N 



In 
x Im 



dx 



dk 



128k 1 / ky/r(k - 2i^x) 2 

o ' xtnn- 1 (k/(2^r)) n V(fc 2 + Ar)(k - Ai^/rx) 
1 



(k - 2i^/rx + (k 2 ~ 6iky/rx ~ 4r(l + x 2 ))uj / '(2y/r(k - 2iy/rx))) 



32iry/T 
N 



In 



Putting the contributions to the last term of Eq. (3.8) in Eqs. ( C3|C4 ) together, we have 



4 f K 4 (p,k,-p,-k) 
3lP(p,r)J fe n(A,r) 



3277^ < 1 



% (A 1-2.008 



(C4) 



(C5) 



Inserting the results in Eqs. (3.5 C2 C5| into Eq. (3.8), we obtain Eq. (3.9) 



Appendix D: Self energies in the Goldstone phase 



We can compute the self-energies associated with the bare connected Green's functions in Eq. (4.7) by a standard 



1/N expansion of the partition function in Eq. (4.5). The self-energies to order 1/N are 
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(Dl) 



Note that all the terms involving 11(0,0) in Eq. (Dl) do indeed cancel with each other. In the last line of each 



expression in Eq. (Dl ) we have extracted the ultraviolet divergences explicitly; the functions F are finite and universal 
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dimensionless functions, which depend only on the ratio pja" 1 . They are given by the following integrals: 
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At small values of p we find 
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whereas for general values of p these functions can be evaluated numerically. Notice that the log(p) terms are present 



only in F\\, and these arise from the integral of the explicitly displayed logarithm in Eq. (D3). 

The imaginary parts of these functions at p = — 16ctq are obtained by analytically continuing the contour integration 
via k — > ke 16 , p — > pe and then taking the limit 9 \ — it: 
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The numerical values of these functions evaluated by Mathematica at a series of values of 9 starting at 
converging to 9 — — 7r, lead to 



(D7) 
0, and 



F aa (-16a%) = 585.93223 + 5 x 10~ 7 i 
F ctA (-16ct 2 ) = -19.0308885 + 3 x 10" 
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-1.5589876- 1.9188981 i, 



(D8) 



in excellent agreement with the values in Eq. (D7|. Note that the imaginary values are restricted to F\\, and it is 



this imaginary contribution which leads to the real part of the pole in Eq. ( 1.5 ). This imaginary value of F\\ is again 



a consequence of the explicitly displayed logarithm in Eq. (D3). Thus the log(p) terms in Eq. (D6| and the oscillatory 



component of the pole are linked to each other, as is also noted in Section [V] 



1 S. Chakravarty, B.I. Halperin, and D.R. Nelson, Phys. Rev. Lett. 60, 1057 (1988). 

2 M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, Phys. Rev. B 40, 546 (1989). 

3 A. V. Chubukov, S. Sachdev, and J. Ye, Phys. Rev. B 49, 11919 (1994). 

4 S. Sachdev, Phys. Rev. B 59, 14054 (1999). 

5 W. Zwerger, Phys. Rev. Lett. 92, 027203 (2004). 

6 T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and M. P. A. Fisher, Science 303, 1490 (2004). 

7 D. Podolsky, A. Auerbach, and D. P. Arovas, Phys. Rev. B 84, 174522 (2011). 

8 R. Coldea, D. A. Tennant, E. M. Wheeler, E. Wawrzynska, D. Prabhakaran, M. Telling, K. Habicht, P. Smeibidl, and 



22 



K. Kiefer, Science 327, 177 (2010). 

S. Sachdev, K. Sengupta, and S. M. Girvin, Phys. Rev. B 66, 075128 (2002). 

J. Simon, W. S. Bakr, R. Ma, M. E. Tai, P. M. Preiss, and M. Greiner, Nature 472, 307 (2011). 

Ch. Ruegg, B. Normand, M. Matsumoto, A. Furrer, D. F. McMorrow, K. W. Kramer, H.-U. Giidel, S. N. Gvasaliya, H. 
Mutka, and M. Boehm, Phys. Rev. Lett. 100, 205701 (2008). 

D. S. Chow, P. Wzietek, D. Fogliatti, B. Alavi, D. J. Tantillo, C. A. Merlic, and S. E. Brown, Phys. Rev. Lett. 81, 3984 
(1998). 

M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, and I. Bloch, Nature 415, 39 (2002). 
I. B. Spielman, W. D. Phillips, and J. V. Porto, Phys. Rev. Lett. 98, 080404 (2007). 
N. Gemelke, X. Zhang, C.-L. Hung, and C. Chin, Nature 460, 995 (2009). 

J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, I. Bloch I, S. Kuhr, Nature 467, 68 (2010). 

W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. I. Gillen, S. Foiling, L. Pollet, M. Greiner, Science 329, 547 (2010). 
M. Endres, M. Cheneau, T. Fukuhara, C. Weitenberg, P. Schaufi, C. Gross, L. Mazza, M. C. Banuls, L. Pollet, I. Bloch, S. 
Kuhr, Science 334, 200 (2011). 

M. Endres, T. Fukuhara, D. Pekker, M. Cheneau, P. Schaufi, C. Gross, E. Demler, S. Kuhr, and I. Bloch, Nature 487, 454 
(2012). 

K. G. Wilson and M. E. Fisher, Phys. Rev. Lett. 28, 240 (1972). 

J. Oitmaa, Y. Kulik, and O. P. Sushkov, Phys. Rev. B 85, 144431 (2012). 

S. Sachdev, Quantum Phase Transitions, Second Edition, Cambridge University Press, Cambridge UK (2011). 
N. H. Lindner and A. Auerbach, Phys. Rev. B 81, 054512 (2010). 

E. Brezin, J. C. Le Guillou, and J. Zinn- Justin, in Phase Transitions and Critical Phenomena, Vol. 6, edited by C. Domb 
and M. S. Green (Academic Press, New York, 1976). 

L. Pollet and N. V. Prokof'ev, Phys. Rev. Lett. 109, 010401 (2012). 

S. Sachdev, in Quantum Theory of Condensed Matter, Proceedings of the 24th Solvay Conference on Physics, Bertrand 

Halperin editor (World Scientific, Singapore, 2010), |arXTv:0901.4103| 

T. Faulkner, Hong Liu, J. McGreevy, and D. Vegh, Phys. Rev. D 83, 125002 (2011). 

F. Denef, S. A. Hartnoll, and S. Sachdev, Phys. Rev. D 80, 126016 (2009). 
S. A. Hartnoll and D. M. Hofman, Phys. Rev. B 81, 155125 (2010). 

R. Abe, Prog. Theor. Phys. 49, 1877 (1973). 

Y. Okabe and M. Oku, Prog. Theor. Phys. 60, 1277 (1978). 



23 



